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The main goal of the KamLAND reactor u e experiment is a search for v £ oscillation using inverse-/? decay reaction 
in 1,000 ton of ultra-pure liquid scintillator. The data-set is 1490.8 days from Mar. 2002 to May 2007. The best-fit 
oscillation parameters are Am^ = 7.58±gj3 (stat.) ± 0.15 (syst.) X 10~ 5 eV 2 and tan 2 6»i 2 = 0.56± ;J° (stat.) l Jg 
(syst.). The statistical significance for reactor P e disappearance is 8.8<r, and an undistorted u e energy spectrum is 
disfavored at > 5a. 
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1. Introduction 
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KamLAND (Kamioka Liquid scintillator Anti-Neutrino Detector) investigates neutrino oscillation parameters 
by observing electron anti-neutrinos (i> e ) from distant nuclear reactors. Previously, KamLAND reported the first 
evidence of D e disappearance [l| and evidence of spectral distortion of the reactor D e energy spectrum Q ■ Furthermore, 
KamLAND published the first investigation of geologically produced anti-neutrinos (geo-neutrinos) from radioactive 
i—,. decay in the Earth [3J. 

In our recent analysis 0, we have enlarged the fiducial volume radius from 5.5m to 6.0m, extended the analysis 
threshold down to the inverse-/? decay energy threshold and reduced the systematic uncertainties in the number of 
target protons and the background. 

,— i '. 2. Detection method and analysis improvements 

> ' 

OO ' 

^ 2.1. Detection method 

v e are detected via the inverse (3 decay reaction, v e + p — > n + e + . The emitted positron annihilates with an 
electron and emits two 0.511MeV 7-rays. This signal appears as a prompt signal and its energy consist of the positron 
kinetic energy and these two 7-rays. The prompt scintillation light from the e + gives a measure of the v e energy, 
_Ep c ~ Eprompt + E n + 0.8MeV, where E n is the average neutron recoil energy, O(lOkeV). The neutron thermalizes 
by elastic scattering in the liquid scintillator and captures on a proton after 207.5 ± 2.8/isec, creating a deuteron 
with 2.2MeV 7 emission. The neutron capture is the delayed signal. By using the delayed coincidence method, the 
backgrounds that do not have timing or space correlations are significantly reduced. 



2.2. Off-axis calibration 

We recently used an off-axis calibration system to investigate the vertex position reconstruction performance. Up 
to 5.5m radius, the vertex reconstruction deviation is within 3cm at any (R, 9, (ft), corresponding to 1.6% fiducial 
volume uncertainty. 12 B/ 12 N spallation products are used to extrapolate this uncertainty up to 6.0m fiducial radius, 
resulting in 1.8% uncertainty. The vertex and energy resolution are also estimated from off-axis and z-axis calibration 
data. The estimated vertex resolution is ~ 12cm/ y/E(MeV) and the energy resolution is 6.5%/ wE(MeV). 



2.3. Likelihood selection 

In Ref.[2|, 5.5m fiducial radius and 2.6MeV prompt energy threshold were applied to the reactor v e analysis. 
To get more statistics and to use the full reactor energy spectrum, a larger fiducial volume and a lower energy 
threshold are required. However, accidental coincidences due to 40 K, 210 Bi and 208 T1 increase near the balloon 
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surface (R = 6.5m) reducing the signal-noise ratio. To maintain high efficiency for v e detection, we apply a likeli- 
hood selection using event characteristics. Firstly we construct a probability density function (PDF) for accidental 
coincidence events, f aC ci(E p , Ed, AR, AT, R p , Rd) by using an off-time (10ms to 20 sec) window. A PDF for the 

*4> fi> e (E p , Ed, AR, AT, R p , Rd) is constructed from a MC simulation. For the E prompt distribution, we select an 

f- 

no-oscillation reactor spectrum including a contribution from geo-neutrinos. A discriminator value, L = -j — '._ 

Jacci-rjv e 

is calculated for all candidate pairs that pass the primary cuts. We choose a selection value L^ ut in E p bins of 0.1 
MeV, where L^ ut is the value which has the maximum figure-of-merit. Table U shows the summary of the systematic 
uncertainties in the analysis. The dominant uncertainties come from the reactor-related ones. 



Table I: Estimated systematic uncertainties relevant for the neutrino oscillation parameters. 





Detector-related(%) 




Reactor-related ( % ) 




Amia 


Energy scale 


1.9 


v e spectra 


0.6 




Fiducial Volume 


1.8 


Ve, spectra 


2.4 


Event rate 


Energy threshold 


1.5 


Reactor power 


2.1 




Efficiency 


0.6 


Fuel composition 


1.0 




Cross section 


0.2 


Long-lived nuclei 


0.3 



2.4. Background 

The dominant background is caused by the 13 C(a,ro) 16 reaction. The primary a source is the decay of 210 Po, 
a decay-daughter of 222 Rn introduced during the liquid scintillator filling. In Ref.Q], a 32% error for the reaction 
to 16 ground state was assigned from uncertainties of 210 Po rate, the cross section to 16 ground state, neutron 
angular distribution and proton quenching. A conservative 100% error was assigned to the excited states since the 
cross section was not well known. To measure the (a,n) prompt spectrum, a detector calibration with a 210 Po 13 C 
source was performed. The expected spectrum was also simulated by the 13 C(a, ?i) 16 cross sections from Ref-f^Q 
and compared with the data. The new cross section measurement Q has a 4% uncertainty in the total cross section, 
but it does not specify the final state separately. Thus the subtraction of the reaction rate of the excited states 
estimated from the cross section in Ref.[6j is needed for this new measurement. The first 16 O excited state emits 
e + e~ , but they cannot escape from the source capsule . In KamLAND, these events are observed at ~ 1.02MeV 
from the annihilation gammas. Therefore, the reaction rate of each state can be measured separately by 210 Po 13 C 
calibration source. The rate of the second excited state was consistent with the calculation jfjj without scaling. For 
the first excited state, the calculated spectrum agreed with the calibration data if scaled by 0.6. After subtraction of 
the excited state cross sections, the ground state was obtained. The reaction rate of the ground state had to be scaled 
by 1.05 to get the calculated ground state to be in good agreement with the data. With these improvements, we 
assign an uncertainty of 11% for the ground state and 20% for the excited states. Finally, the number of background 
events due to 13 C(a,n) 16 is estimated to be 182 ± 21.7 events. Tabic HT1 shows the summary of the backgrounds. 



3. Results 

We observe 1609 events in our data-set. With no v e disappearance, we expect 2179 ± 89 events from the reactors. 
The background in this analysis is 276.1 ± 23.5 events excluding the geo-neutrino signal. Figure [1] shows the prompt 
energy spectrum of selected v e events and the fitted backgrounds. The unbinned data is assessed with a maximum 
likelihood fit to two-flavor neutrino oscillation, simultaneously fitting the geo-neutrino contribution. The method 
uses the absolute time of the event and accounts for time variations in the reactor flux. Earth matter oscillation 
effects are also included. 
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Table II: Summary of estimated backgrounds 



Background 


contribution 


Accidentals 


80.5 


± 


0.1 


8 He/ 9 Li 


13.6 


± 


1.0 


Fast neutrons and Atmospheric v 


< 


9. 





13 C(q, n) 16 G.S. 


157.2 


± 


17.3 


13 C(a, n) ie O 12 C(n, n 7 ) 12 C (4.4 MeV 7) 


6.1 


± 


0.7 


13 C(a, n) le O 1st excited state (6.05 MeV e + e") 


15.2 


± 


3.5 


13 C(a, n) ie O 2nd excited state (6.13 MeV 7) 


3.5 


± 


0.2 


Total 


276.1 ± 23.5 




E p (MeV) 



Figure 1: Prompt energy spectrum of Ve. candidate 
events. All histograms corresponding to reactor spec- 
tra and expected backgrounds incorporate the energy- 
dependent selection efficiency (top). The shaded back- 
ground and geo-neutrino histograms are cumulative. 
Statistical uncertainties are shown for the data. The 
band on the blue histogram indicates the event rate sys- 
tematic uncertainty. 



The best-fit is shown in Figure [T] The best-fit oscillation parameters are Am^ = 7.58^ (stat.) ± 0.15 (syst.) 
x 10~ 5 cV 2 and tan 2 #i 2 = 0.56+q;q° (stat.) ig]o6 (syst.). A scaled reactor spectrum with no distortion from 
neutrino oscillation is excluded at more than 5 a. Figure [2] shows the allowed contours in the neutrino oscillation 
parameter space. Only the so-called LMA-I region remains and the LMA-II and LMA-O regions are disfavored at 
> 4(7. A combined oscillation analysis of KamLAND and other solar neutrino experiments under the assumption of 
CPT invariance gives Am^ = 7.59 ± 0.21 x 10~ 5 cV 2 and tan 2 6>i 2 = 0.47±o.o5- Tne Am 2i parameter is strongly 
determined by the KamLAND experiment. 

The Lo / E distribution is suitable for illustrating the oscillatory behavior, see Figure[3] In this figure, the data 
and the best-fit spectrum are divided by the expected spectrum without oscillation. L is the flux-weighted effective 
baseline and is 180 km. Two oscillation periods are seen, the KamLAND data is in excellent agreement with the 
oscillation scenario. 

To determine the number of geo-neutrinos, we fit the normalization of the D e energy spectrum from U and Th 
decay chains simultaneously with the neutrino oscillation parameters using the KamLAND and solar data. Fixing 
the Th/U mass ratio from planetary data to 3.9, the best-fit value is 73 ± 27 events. 



4. Conclusion 

The KamLAND experiment has determined a precise value for the neutrino oscillation parameter Am 21 and 
stringent constraints on #i 2 using reactor anti-neutrinos from 55 Japanese reactors over a five years period. The 
observed neutrino spectrum is distorted and the significance of the spectral distortion is > 5<r and no-spectral 
distortion is strongly disfavored. A two-neutrino oscillation analysis of the KamLAND data with 0.9 MeV threshold 
gives tan 2 6 12 = 0.56±g;J° (stat.) t° m (syst.) and Am^ = 7.58 (stat-) ±0.15(syst.) x 10" 5 eV 2 . The previous 
result [3] was tan 2 #i 2 = 0.46iQ 2 g and Am^ = 7.9 ^3 x 10~ 5 e V 2 - The uncertainties of neutrino oscillation 
parameters are significantly reduced and the systematic and statistical uncertainty are now comparable. The LMA- 
II and LMA-0 which were both allowed in Ref. [2j are now disfavored at > 4 a and only LMA-I remains. These 
results are also published in Ref. Q . 
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Figure 2: Allowed region for neutrino oscillation pa- 
rameters from KamLAND and solar neutrino experi- 
ments. The side panels show the Ax 2 -profiles for Kam- 
LAND (dashed), solar experiments(dotted) and the com- 
bination of the two(solid). 
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Figure 3: Ratio of the background subtracted v £ spec- 
trum to the expectation without oscillation as a func- 
tion of Lo /E. Lo is the effective baseline taken as a flux 
weighted average (Lo = 180 km). The energy bins are 
equal probability bins of the best-fit including all back- 
grounds. The blue curves are the expectation account- 
ing for the distance to the individual reactors, time de- 
pendent flux variations and efficiencies. The error bar 
shows the statistical uncertainty only. 
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